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Abstract. 

The boundary of the Siegel disk of a quadratic polynomial with an irrationally 
indifferent fixed point with the golden mean rotation number has been observed to 
be self-similar. The geometry of this self-similarity is universal for a large class 
of holomorphic maps. A renormalization explanation of this universality has been 
proposed in the literature. However, one of the ingredients of this explanation, the 
hyperbolicity of renormalization, has not been proved yet. 

The present work considers a cylinder renormalization - a novel type of 
renormalization for holomorphic maps with a Siegel disk which is better suited for 
a hyperbolicity proof. A key element of a cylinder renormalization of a holomorphic 
map is a conformal isomorphism of a dynamical quotient of a subset of C to a bi- infinite 
cylinder C/Z. A construction of this conformal isomorphism is an implicit procedure 
which can be performed using the Measurable Riemann Mapping Theorem. 

We present a constructive proof of the Measurable Riemann Mapping Theorem, 
and obtain rigorous bounds on a numerical approximation of the desired conformal 
isomorphism. Such control of the uniformizing conformal coordinate is of key 
importance for a rigorous computer-assisted study of cylinder renormalization. 
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1. Introduction 

Over the years, one particular mathematical tool - renormalization - has become very 
instrumental in problems involving universality in a class of maps or flows. 

It is possible to give a very concise definition of renormalization in the dynamical 
setting: A renormalization of a dynamical system X acting on a topological space T is 
a rescaled first return map of a subset of T. 

To prove universality, one usually introduces a renormalization operator on a certain 
functional space, and proceeds to show that a) this operator has a fixed point, and b) 
the operator is hyperbolic at this fixed point. 

This approach has led to dramatic successes in the rigorous explanation of the 
universality in several areas of complex and real dynamics, most notably, universality of 
unimodal maps through the works of Sullivan (1986, 1992), McMullen (1998, 2000), 
Lyubich (1995, 1999) and others, universality of critical circle maps: de Faria 
(1992, 1999), de Faria and de Melo (1999, 2000), Yampolsky (2001, 2002a, 2002b) and 
others, and universality of "critical" Hamiltonian flows: Koch (2002, 2004a, 2004b) and 
others. 

However, the renormalization picture for one of the most interesting objects in 
one-dimensional complex dynamics - Siegel disks - is not complete yet. Recall that a 
Siegel disk is the maximal linearization domain of a holomorphic map with a linearizable 
irrationally indifferent fixed point. 

A classical theorem of Siegel implies that for the quadratic polynomial 



linearizable. The Siegel disk Ag* of fg* is a quasidisk whose boundary contains the 
critical point: 1 G dAg* (cf e.g. (Douady 1987)). It has been observed numerically 
by Manton and Nauenberg (1983) that the boundary of the Siegel disk of fg* is 
asymptotically self-similar. They found that the Fibonacci iterates of the critical point, 
fg*(l), n = 1,2, 3, 5, 8, . . ., approach the critical point asymptotically along two straight 
lines separated by an angle 2a 107.3 degrees. At the same time, the distance between 
successive iterates decreases at a geometric rate: 



fa = e 2ni9 *z(l - 0.5z) 



(1.1) 



with 9* = (v 7 ^ - l)/2 (the inverse golden 




lim 



(1.2) 



/*n(l)_l 



with A ~ —0.7419. Thus, the geometry on the boundary of the Siegel disk repeats itself 
around the critical point after a reflection with respect to a certain line passing through 
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point 1, and a rescaling. The self-similar geometry is transported by the dynamics to 
other points on dAg* . 

Manton and Nauenburg also conjectured that the parameters a and A, as well, 
as the boundary curve of the Siegel disk at the critical point and its preimages, are 
universal for a large class of maps with a fixed point of multiplier e 2md * at the origin. 
Widom (1983) considered a renormalization operator defined on pairs of such commuting 
holomorphic maps, and conjectured the existence and hyperbolicity of a fixed point in 
a space of such pairs. He also obtained numerical evidence that such a fixed point does 
exist. 

In the 90 's there were several results that used renormalization to explain this self- 
similar structure. Stirnemann (1994) was able to construct a computer-assisted proof 
of the existence of a renormalization fixed point for holomorphic maps with a golden 
mean Siegel disk. More recently, McMullen (1998) proved the asymptotic self-similarity 
of the Siegel disk of the golden mean quadratic polynomial fg*. He has demonstrated 
that successive renormalizations of the golden mean quadratic converge to a fixed point. 
However, the renormalization for Siegel disks based on commuting pairs has not yet led 
to a proof of the hyperbolicity of the renormalization operator. 

The main difficulty with an attempt to prove hyperbolicity through the 
renormalization of commuting pairs is that the space of commuting pairs does not 
possess the natural structure of a Banach manifold in which renormalization is an 
analytic operator. Such pairs constitute a topological space whose properties are far 
from simple. 

Yampolsky (2002a) introduced a cylinder renormalization of analytic maps with a 
fixed point. The novelty of this renormalization is that it is formulated in terms of an 
analytic operator on a certain Banach space of analytic functions. Such an operator is 
more suitable for a proof of hyperbolicity. The present work fills out several steps in a 
computer-assisted implementation of such a proof. 

Cylinder renormalization has been adopted for holomorphic maps with a Siegel disk 
in (Yampolsky 2005). We will continue with a brief summary of Yampolsky's procedure. 
Define A p to be the Banach space of all bounded analytic functions on the open disk 
of radius p > 1 around zero, D p C C, equipped with the sup norm. C p will denote the 
subspace of A p which consists of all / G A p for which is a fixed point, /(0) = 0, and 
1 is a critical point, /'(l) = 0. 

Given an / G C p , suppose that for some n G N there exists a simple arc / connecting 
a fixed point a n of f n with the fixed point such that f~ n (l) is again a simple arc 
intersecting I only at points a n and 0. We will call the interior of the region bounded by 
the curves / and f~ n (l) the fundamental crescent Cp of cylinder renormalization if the 
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Figure 1. Schematics of the cylinder renormalization. 



iterate f~ n \c fn is univalent and the quotient of Cfn U f~ n (Cfn) \ {0,a n } by the iterate 
f n is conformally isomorphic to C/Z. 

Furthermore, let us assume that Rfn, the first return map of Cfn, has a critical 
point c G Cfn. The conformal isomorphism of Cfn and C/Z, normalized so that it maps 
c to 1, will be denoted by $. The map n R cy i{f) = e o $ o Rfn o o e _1 (see Fig. 
1), where e(z) = e~ 2mz , will be called a cylinder renormalization of / with period n if 
nRcyi(f) £ C p / for some p' > 1. Naturally, this definition should be independent of the 
particular choice of the fundamental crescent. In fact, any other fundamental crescent 
C'f„ with the same endpoints as Cfn, and such that C/« U Cj n is a topological disk, 
produces the same n R C yi{f) (cf (Yampolsky 2002a)). 

It is clear from this definition that the knowledge of the conformal isomorphism of 
Cfn and C/Z is of a key importance for building a cylinder renormalization of a map. In 
the remaining part of the paper we will describe how to construct such an isomorphism 
for an analytic map of degree d > 2 with an irrationally indifferent fixed point at 0: 



/(0) = and f(0) = e 2md . 

We will further assume that 6 is of bounded type, that is it admits a continued 
fraction expansion 

° = r T 1 i » (1-3) 

(which we will abbreviate as = [r , ri, r 2 , . . .]) with 

suprj < oo. 

The rational number corresponding to the n-th truncation of this expansion is 
commonly denoted by q n /p n : 

— = [r ,ri,...,r n ]. 

Pn 
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By the classical result of Siegel (1942), for such an / there exists a maximal open 
connected neighborhood U of 0, called the Siegel disk, on which the action of / is 
conformally conjugate to the rigid rotation z i— > e 2 z of the open unit disk. 

We will now focus on the construction of a cylinder renormalization for analytic 
maps with a Siegel disk. X. Buff (unpublished) has suggested a construction of a 
fundamental domain for a quadratic polynomial as the interior of a region bounded 
by internal and external rays, logarithmic spirals and equipotentials. However, the 
exponential density of the Julia set makes a computer implementation of some of these 
curves practically difficult. 

Here, we will make a different choice of a fundamental crescent for functions with 
a Siegel disk. The curve l n will consist of two parabolas: one passing through points 
and /?»+2+9n(l) ) the other - through point /9™+2+9«(l) an d a repelling fixed point 
a qn . Such two parabolas are uniquely defined after one specifies their common tangent 
line at point fi n + 2 +i n (l^. The slope of this line can be chosen in a convenient way. Of 
course our, choice of the boundary l n is rather arbitrary, but has the virtue of having a 
simple analytic form. This simple choice of boundaries can be used for a general analytic 
function in C p , and not only for quadratic polynomials. 

We shall take the interior of the region bounded by the curves l n and f~ 9n (l n ) as 
a candidate for a fundamental crescent Cfq n (here we chose the inverse branch of f~ Qn 
that maps into itself). One has to verify that the simple curve f~ q "(l n ) intersects l n 
only at the endpoints. We will momentarily defer the explanation of this verification. 

As we have already explained, the key step in the cylinder renormalization of a 
function / is a construction of a properly normalized conformal isomorphism $. It 
should be noted however that this uniformizing change of coordinates can not be written 
down explicitly, and, as we will show, requires a solution of a Beltrami equation which, 
in turn, can not be found "by hand" . One, therefore, may try to find $ numerically and 
provide bounds in this numerical solution: Below, we will demonstrate how this can be 
done for the golden-mean quadratic polynomial fg*, and n — 3. 

2. Uniformization of a fundamental crescent of cylinder renormalization 

In this Section we will outline the steps in the construction of the conformal isomorphism 
$ (see Fig. |2] for the schematics of our procedure). Most of the details of this 
construction are independent of the choice of /, however at the end of this Section 
we will specify our fundamental crescent for the golden-mean quadratic polynomial and 
n — 3. 
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Figure 2. Mapping g 1 1 o r 1 of the interior of the fundamental crescent C/ 
the punctured plane with a slit C* \ IR + . 



onto 



Introduce for all z G Cfq n a new coordinate £ through 

r(0 = t 



5n 

p ia^-\-b 



[2.1) 



where normalizing constants a and b are chosen so that |r _1 (/ 9n+2+,?n (l))| = 1 and 
T _1 (/ 9n+2 (l)) = 0. The choice of of this coordinate is motivated by the fact that r _1 
maps the interior of the fundamental crescent conformally onto the interior of an infinite 
vertical closed strip S, whose width is comparable to one. This coordinate has been 
suggested in (Shishikura 1942) in a similar context. 

Next, following (Shishikura 1942), define a function g n from U == {u + iv G C : 
< Rew < 1} to S by setting 

g n (u + iv) = (1 - u)T-\f- qn ( ln {v))) + UT-\ ln {v)) 

(below, we will specify the choice of parametrization ) n : I h l n ). Let do be the 
standard conformal structure on C, and let a = g^cr be its pull-back on U. Extend 
this conformal structure to C through a = (T k )*o~ on T~ k (U), T(w) — w + 1. 
By the Measurable Riemann Mapping theorem (see Section 0J) there exists a unique 
quasiconformal mapping g : C i— > C such that g*o~Q = a, normalized so that ^(0) = 
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and g(l) = 1. Notice that g o T o g 1 preserves the standard conformal structure: 
(goTog- 1 )* a = (g' 1 )* o T* og*a = (f'foTV = {g^fa = (<f)~V = a 



therefore it is a conformal automorphism of C, i.e. an affine map, and in fact, a 
translation, since T does not have fixed points in C. Moreover g o T o <7 _1 (0) = 1, 
Therefore g o T o g^ 1 = T. 

By the definition of g n , g~ x o r _1 o f qn o r = T o g n _1 on the image of Z n by r _1 . 
Therefore, the map <p = g ofi^ 1 , defined on 5, can be continuously extended to all of C. 
Since is conformal on the interior of S, § <p(w)dw = over any closed contour 7 in C, 
and by Morera's theorem is analytic on all of C. Clearly, $ = o r _1 is the desired 
conformal isomorphism of C/?n and C/Z. 

Next, define g = e o g o e _1 . Since 

9s(e(w)) = e(w) <fo(w) 
^(eM) e (iu) fiWw)' 

the 1-periodic function <? is a solution of the Beltrami equation 

gw = fj>g w , p. = {g n )w/ {g n ) w 

whenever g is a solution of 

9s = t*9z, K z ) = (z/^He' 1 ^))- (2-2) 

Thus, we have reduced the problem of finding eo$ = joeo g^ 1 o r _1 to that of 
finding the properly normalized solution of the Beltrami equation 

/ x z {9n)w{e~ 1 {z)) . , 

9 ' = " 9 " " w = iao.(«-'W) ( ] 

on the punctured plane C*. The issue of existence of a solution of (|2.3|) is addressed 
by the famous Measurable Riemann Mapping theorem, also known as Ahlfors-Bers- 
Boyarskii theorem (see (Ahlfors and Bers 1960) and (Boyarskii 1957)): 

Theorem 1. (Ahlfors and Bers, Boyarskii). Let fi G Loo(C) satisfy WfiW^ < K < 1. 
Then there exists a unique map g^ : C — > C such that g%{z) / 'g£(z) = ^{z), g^ fixes 0, 1 
and 00, and g^ is a (1 + K)/{1 — K)-quasiconformal map. 

Recall, that there are two equivalent definitions of a quasiconformal map, the 
"geometric" and the "analytic" one. Here we will give only the "analytic" definition 
which is shorter: an interested reader is referred to (Ahlfors 1966), (Bers 1977), 
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(Markovic 2003) and (Boyarskii and Iwanier 1974) for a fuller account of the basics 
of the theory. 

Given an open set Q C C, a map / : Q — > f(Q) C C is said to be absolutely 
continuous on lines (ACL) in a rectangle R C Q with sides parallel to the x and y axes 
if / is absolutely continuous on almost every horizontal and vertical line in R. The map 
/ is ACL on Q, if it is ACL on every rectangle in Q. Partial derivatives f z = (f x — if y )/2 
and fz = (fx + ify)/2 of such map exist a.e. in fl. 

Definition 2. A homeomorphism f : Q — > f{Q) is K-quasiconformal if and only if the 
following holds: 

i) f is ACL on every rectangle in Q, 

a) \M < fril^l «- e - in ^- 

The complex dilatation of a quasiconformal homeomorphism / is fif(z) = 
fz(z)/ f z (z). A mapping / with a prescribed complex dilatation \i solves the following 
Beltrami equation: 

Mz) = (ji(z)f z (z), (2.4) 
and its existence is the subject of Theorem ^ 

3. Statement of main results 

Most of the present work deals with a constructive implementation of the original proof 
of the Measurable Riemann Mapping Theorem due to L. Ahlfors and L. Bers. One of 
our main objectives will be to present rigorous bounds on how far a numerical solution 
of ()2.3|) lies from the actual one. We will give a detailed derivation of the necessary 
bounds in Sections loTfTUl At the moment, however, we will proceed to summarize our 
results. 

Theorem 3. Let /i G Loo(C) and an integer p > 2 be such that WfiW^ < K < 1 and 
KC P < 1, where 

C p = cot 2 (7r/2p). 

Assume that \l = v + 77 + 7, where v and r\ are compactly supported in and j(z) is 
supported in C\Dr. Furthermore, let 77 be in L p (Dr) and \\rj\\ < 5 for some sufficiently 
small 5. Also, let h* e L P (C) and e be such that B p (h*,e), the ball of radius e around 
h* in L P (C), contains B p (T v [h*], C p e') , with 

e = 5 ess sup \h* + 1| + Ke. 



Renormalization for Siegel disks and Measurable Riemann Mapping Theorem 



8 



Then the solution g 1 * of the Beltrami equation g% = fj,g% admits the following bound: 

C[\g»(z)\ + Ae'\z\^] 1+2/p 

R 

where g*(z) = P u [h*](z) + z and 



"(*) - 9 :(z)\ < Ae'\z\^ + U^i^l^ (3.1) 

it;4/ P -C[|^(z)|+Ae'^|i- 2 /P] 2 ^' 



A = 1 



7T 



VP 



3 2< ^ 2 25 f-l 



"29(2 -g) 36 g-1 



C = * 1 '*Ak !?££ ™ p } , r = inf • 
r 2 /p(l - XCp) { 2 ee| 2 |=iJ} u W|J 

We have used the bounds from Theorem |3] to uniformize a fundamental domain 
for the cylinder renormalization of period 3 of the golden mean quadratic polynomial. 
Before we state our findings, we would first like to specify our choice of a fundamental 
crescent Cp. 

It is convenient to parametrize the boundaries l n and f~ qn (l n ) of Cfq n by the radial 
coordinate in C. Of course, this parametrization is far from unique. Here, we will 
mention our choice for fg*(z) = e 27Tte * z(l — 0.52;) and n — 3, motivated by our numerical 
experiments: 

A 3 (r) = ( ( x ( r )'^ x ( r ) 2 + Bx ( r )), r<f, 

\(Cy(r) 2 + Dy(r) + E,y(r)), r > f, 



where 



Rejff(l) 

x(r) = W T(r)l 

v(r) _ i m f n (1) \as-nm\ + \nm-nr) . T^-mm 
\^-nm + \nm 



T(r) 



and f is defined through T(f) = |/g«(l)|. Constants A, B, C, D and E are fixed by 
the conditions 0, fg}(l),a 3 G Z3, together with the requirement that the slope of the 
common tangent line to both parabolas at point fg*{l) is equal to —1.7. 
Define the following function on C*: 

Sn(r, 0) = U-<f>) + r- l (f-^(X n (r))) + fl - ^(-0) - |-) r~\K{r)\ 

where — tt < < ir, and 77 is the Heaviside step function (we have adopted the convention 
?7(0) = 1). It is clear that g n o e _1 (re^) = g n (r, 0). Therefore if fi is as in (|2.3jl . then 

„( re <«) = c 2i<i> rd r9n(r, 4>) +id 4 ,g n (r, 0) 
rd r g n (r, 0) - id^g n (r, 0) 
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Figure 3. The fundamental domain for the cylinder renormalization of the quadratic 
polynomial f ,{z) = e 2m9 * z(l - 0.5z), 9* = (\/5- l)/2, of period 3 (a) and the Siegel 
disk of the cylinder renormalization of fg* of period 3 (b). 



on C*. 

One also needs to verify the empty intersection of the boundaries of the fundamental 
domain. To this end, choose a sufficiently large R and a sufficiently small g, and consider 
the arc X n ([g, R]). Cover this arc by a finite collection {Bi} i=1 of open balls in C (see 
Appendix for a discussion of evaluation of analytic functions on open disks, referred to 
as "standard sets in C) and verify that f~ qn (Bk)U{Bi}f =1 — 0, k — 1 ... N. Verification 
of the empty intersection of the boundary curves for small, r < g, and large, r > R, 
values of the parameter is not needed: Rather than using parabolas, we will complete 
the boundaries of the fundamental crescent for these parameter values with pieces of 
rays and logarithmic spirals together with their preimages (see Section|U]for details). We 
will refer to this construction simply as a "completion" of the fundamental crescent. We 
have performed the verification of the empty intersection for the golden mean quadratic 
polynomial fe*(z) = e 2md * z{l — 0.5z), and n = 3. Although most of the results of 
this paper are valid for any analytic / G C p , our Main Example, that serves purely 
demonstrative purposes, will be stated only for these particular data. 

Main Example. Let f e *{z) = e 2md * z{l - Q.hz), where 9* = (V5 - l)/2, and let A 3 be 
as in (jj.jjj) . Then, there exists a completion of the fundamental crescent bounded 
by the arcs X 3 ([g,R}) and fg^(X 3 ([g,R})) with g = 3.1 x 10~ 3 and R = 9.8 x 10 3 ; and 
a finite Fourier series whose dilatation is compactly supported in B>r such that the 
solution of the Beltrami equation (jg.,^ on C satisfies 

I MM <M<A>\ ii-2/p , C Pg»(*)l + Ae'\zf- 2/p ] l+2/P 
\g^(z) — gJz) \< Ae \z\ tP H ttt-, 

R^ - C [\g,(z)\ + Ae>\z\^/pf p 
forp= 2.1 and A < 7.11, e' < 2.66 and C < 6.14. 
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All bounds on the constants reported in the preceding Example have been found 
with the help of a computer (see the following Sections for a detailed discussion of 
the ingredients of the computer-assisted proof). To this end, we have written a set of 
routines in the programming language Ada 95 (cf (Taft and Duff 1995) for the language 
standard). We have parallelized our programs and compiled them with the public version 
3.15p of the GNAT compiler (Ada Core Technologies 2004). The programs (Programs) 
have been run on the computational cluster of 64 2.2 GHz AMD Opteron processors 
located at the University of Texas at Austin. 

We believe that out constructive bounds on the solution of the Beltrami equation 
arising in cylinder renormalization are a major step towards rigorous computer-assisted 
study of this operator. Apart from that, a constructive proof of the Measurable Riemann 
Mapping Theorem is interesting in its own right, and should find other applications in 
one-dimensional complex dynamics and geometry. 

4. A constructive proof of the Measurable Riemann Mapping Theorem 

The classical proof of Theorem ^ uses two integral operators in an essential way; the 
first operator being the Hilbert transform 

T[h](z) = ^-lini J ! J%- 2 d£Ad£, 
the second — the Cauchy transform 

™M = s// c <&*** 

Some of the properties of the Hilbert transform are described in the celebrated 
Calderon-Zygmund lemma (see (Calderon and Zygmund 1956)): 

Lemma 4. (Calderon and Zygmund). The Hilbert transform is a well-defined bounded 
operator on L p (C) for all 2 < p < oo: for every such p there exists a constant C v such 
that lirf/iJH < C p ||/i|| for any h G L p (C), and C p — > 1 as p — > 2. 

For compactly supported differentials one can state an extended version of the 
Ahlfors-Bers-Boyarskii theorem: 

Theorem 5. Let /j, G Loo(C) be compactly supported and satisfy WfJ-W^ < K < 1. Then 
for every p > 2 such that C P K < 1, the operator h — > T[fi(h + 1)] is a contraction on 
L P (C) with a unique fixed point h* . Moreover, the solution of the Beltrami equation, g^, 
is given by 

g" = P[fi(h* + l)]+id, (4.1) 
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and is such that g(0) = 0, g is continuous, has distributional derivatives, and g z — 1 G 



In the rest of the paper the operators T\ji(- + 1)] and P[fi(- + 1)] will be abbreviated 
as T M and P M , respectively. 

The above version of the Ahlfors-Bers-Boyarskii theorem for compactly supported 
differentials provides the building blocks for the standard proof of Theorem ^ (see 
(Ahlfors 1966)). Here we will present a brief sketch of this proof. 

Let \i = v + 7 where 7(2) = whenever z G Dr, the disk of the radius R around 
0, and v(z) = whenever z G C \ Or. Let g v be the solution of the Beltrami equation 
g\ = vg v z . First, define j3 by setting 



which is compactly supported in a neighborhood of zero, and let g" solve the Beltrami 
equation for (3. Then, g p (z) = R 2 /g' p {R 2 /z) solves g% = Pg%. Finally, = g? o g v 
solves g$ = fj,g%. 

The Beltrami equation arises naturally in problems of uniformization of domains 
in C. Many uniformization problems require rigorous estimates on the solution of 
the Beltrami equation. In some cases one can obtain rough, but sufficient estimates 
"by hand" (cf (Shishikura 1942) for an example similar to one considered here); 
often, however, one requires much tighter bounds on the solution. Although there 
are efficient numeric algorithms for a Beltrami equation with a compactly supported 
/i (cf (Daripa 1992), (Daripa 1993), (Daripa and Mashat 1998) and (Gaydashev and 
Khmelev 2005)), rigorous estimates on how far such numeric solutions are from a true 
solution have not been yet provided in literature. Our main objective will be to describe 
how one can obtain such estimates on an approximate numerical solution. In doing 
this, we will consider a Beltrami equation for the differential (|3.3|) that arises in the 
abovementioned problem of cylinder renormalization. 

5. The Calderon-Zygmund constant 



L P (C). 



H i) — ! - v - 

9z 

Next, one obtains the solution g@ of g§ = /3g%: Define 




(4.2) 



P(z) = (3(R 2 /z)z 2 /z 2 , 



(4.3) 



In our following exposition we will require an estimate on the Calderon-Zygmund 
constant. 
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Lemma 6. The Calderon-Zygmund constant C p satisfies 

C p < cot 2 (7r/2p). 

Proof. For f(x) G Cq(R) define the one-dimensional Hilbert transform 

f(y) 



H[f](x) = P.V.- r J^-dy. 

7T y-x 



-ooV-% 

Application of the Calderon-Zygmund lemma to this singular transform yields 



\H[f\\l<A P 



(5.1) 



As it has been shown in (Grafakos 1997), the best bound on A p for 2 < p < oo is 
cot(7r/2j9). Next, for g G C%(C) define 



T*\q](z 



i 

Air 




By a classical result (cf (Ahlfors 1966)) T[g] = -T*[T*[g}} for g G C 2 (C), and this 
equality can be extended to L P (C). Moreover, T*[g] G L P (C) whenever g G L P (C), and 



T*\q](z 



1 f w fl f°° g(z + re**) - g(z - re i(t> ) 



Hence, 



\T*[g]\\ p <-max 



7i 



dr I e 



g(z + re 1 



g\z — re 



ii>\ 



-dr 



o 



max 



max 



71 



g(z + re 1 *) — g(z — re 1 *) 



dr 



o 



g(z + re 



i<f>\ 



-dr 



max 



l^lllp: 



(5.2) 



where g^z) = g(ze 1 *). The L p -norm in ()5.2)1 is two-dimensional, unlike in (jQJ) . 
However, since 



\H 




\H\g^}{z)\ p dxdy 



dy I \ H [g<p]( x + iy)\ p dx 

"2tt poo 



< / dyA* / \ 94> {x + iy)\ p dx = Al 



o jo 



\g(pe^ + *)\ P pdpdi) 



POO 



\g(pe^)\ P pdpd^ = Al\\gf p 



(5.3) 
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one can get an estimate on the operator T*[ ■ ] in terms of the constant A p from the 
one-dimensional bound (|5.1j) : 

||T*^]|| p <A p y| p <cot(7r/2p) \\g\\ p . 

Finally, 

i™i| p <cot 2 (7r/2p) \\g\\ p . 

□ 



6. Estimates on the fixed point of the operator T M 

Theorem El suggests the following approach to the rigorous computer bounds on the 
solution of the Beltrami equation on the Riemann sphere: Write /i as \l = k + 7 with k 
and 7 supported in the neighborhoods of zero and infinity, respectively, find numerically 
an approximate fixed point of the operator T K for a n with a large compact support, 
introduce "standard sets" in L p (C) (we will postpone the discussion of standard sets for 
later, right now one can think of a standard set in L P (C) as a ball in L P (C)), choose a 
standard set that contains our numerical approximation, and show that this standard 
set contains its image under the operator T K . Then the contraction mapping principle 
guarantees that the same standard set contains the fixed point h* of the operator T K . 
Finally, compute P K [h*} + id, find bounds on this operation, and estimate the error due 
to ignoring the differential 7. 

We will start with a Beltrami differential fi G L £X) (C), ||/i||oo < K < 1, given by 
fi = k + 7 where k and 7 are compactly supported in a large disk B# and a small 
neighborhood of infinity C \ Dr, respectively. We will further assume that k consists of 
a "known" part v (for example, representable as a finite Fourier series, see Appendix), 
and an "error" 77 with a bounded L p norm: k = v + rj, \\f]\\ p < 5 for some sufficiently 
small 5. 

Let h* G L p (C) be an approximate fixed point of the operator T u on L P (C). Also, 
let B p (h*,e) denote a ball around h* in L p (C) of radius e, i.e the set of all functions 
h — h* + f in L p (C) such that ||/|| < e. By contraction mapping principle the ball 
B p (h*,e) contains a fixed point of T u+V if T u+v [B p (h*, e)] C B p (h*,e). The following 
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estimate is straightforward: 

II WO - Uh*)\\ p = \\T[ V (h*+ !) + (!/ + V)f]\\ p < C p \\v(h*+ l) + (y + rj)f\\ t 



<c n 



(ess sup \h*+ l\ p ) p WvWr, + ( ess SU P \ v + ^1 



\r}\\ „sup \h* + 1| + (ess sup |z/+77| p )p 



(6.1) 



where we have used the Holder inequality and the Calderon-Zygmund lemma. One can 
use the fact that for a fixed p < oo and any measurable set E with a finite measure 



ess sup \f\ p = lim 



l/l 



pn 



1/n 



lim 

n— >oo 



P 

E,pn 



ess sup | / 1 

E 



to rewrite (|6.1|) as 



\T V+V (h) - T v (h*)\\ =C P 



5 sup |/i* + 1| + Ke, 



(6.2) 



The above estimate is just one of the several possible ways to bound the difference 
T u+V (h) — T u (h*) in terms of smallness of rj and /. The choice of norms that measure 
this smallness, in particular ||^|| p (as opposed to a seemingly more natural H^Hoo, is 
specifically adopted to our model uniformization problem (see Section^). 

The estimate ()6.2|) gives a bound on the image of a standard set under T v+r] : If 
one can show that the ball B p (T u (h*), C p e'), e' — S sup D \h* + 1| + Ke, is contained in 
B p (h*, e) for some p such that KC P < 1, then the fixed point of the operator T v+rj lies 
in B p (h*, e). 



7. Estimates on the solution of the Beltrami equation on a compact 



Define 

g:Xz) = z + P v [h%z). 

The proximity of the exact solution of the Beltrami equation on the compact to 
g%(z) is described in the lemma below. 

Lemma 7. Let the Beltrami differential v +rj be compactly supported in and satisfy 
\\u H-^lloo < K < 1. Also, given an integer p > 2, such that KC P < 1, let \\rj\\ < 5. 
Furthermore, assume that the fixed point h of the operator T u+V lies in an ^-neighborhood 
ofh* G L P (C), i.e. he B p (h*,e). Then 



\g^(z)-g:(z)\<A 



5 sup \h* + 1| + eK 



il-2/p 



(7.1) 
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A 



71 



1/p 



4:r 4 r + f£ 3 2* - 23 2<? " 2 + 



(I) 



15 



(7.2) 



29(2 -q) 36 " q-l 
and q is the conjugate exponent of p. 

Proof. Recall that the solution of the Beltrami equation is given by formula ()4.1|) . A 
rather standard method (see, for instance (Carleson and Gamelin 1991) ) enables ones 
to get easy Holder estimates on \g u+r '(z) - g"(z)\ = \P V+V [h](z) - P v [h*\{z)\. First, 

\g v+n {z) - £(z)\ = \P[v(h* + l) + (v + T,)(h- h*)](z)\ 
1 



2tt 




iv(0(h*(0 + 1) + HO + q(fl)/(fl]* ^ A d( _ 



<- Uh* + l) + (v + T,)(h -h*)\\ p 
71 1 



\z\ q rdrd<p 



1/9 



where q is the conjugate exponent of p and / = h — h*. Next, split the last integral 
into three pieces: D\ = {£ 6 C : |£ — z| < |^|/a}, .D 2 = {£ G C : |£| < |z|/a}, for some 
a > 2, and _D 3 = D/j \ -Di U Z^2- The first integral over D\, where |£| > (a — 1) |^|/a, can 
be estimated in the following way: 



rdrd<j) < 



a — 1 



1 



-rdrdcf) < 



a 2 *- 2 2tt 
(a - 1)9 2 - g' 



1 2-9 



The integral over _D 2 is controlled in a similar way. Next, let D f 3 = {£ G C : 
|£ - «/2| < &|z|} \AU D 2 , for some 6 > 1/2 + 1/a, and £>£ = {£ G C : |£ - */2| > 
6|2;|, |£ — z/2\ < R + \z\/2}, then the integral over D 3 can be bounded as follows: 





z 


[Jd 3 





rdrd(j) < 



|z| 9 r<ir<i0 



+ 



< 



D 



< a 2q 7x [ b z — ^ ) \z 



\z\ 


7 rdrd(f) 


1 z 
\ a 




2 

a 


1 



+ 27T 



3 2 



(r 2 _ B_ )ff 



|2-g 



71 



(b 2 ~ j) 
q-1 



IN 1-9 



1 2-9 



Therefore, 

\P[rj(h* + l) + (u + r])(h - h*)](z)\ < A \\ V (h* + 1) + (v + r})(h - h* 
where 



.4 



7T 



1/p 



7 2g-2 



(a - 1)«2- q 



+ a 2 % 2 - 2a 2q - 2 + 



(b 2 ~ j) 
9-1 



1\W 
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The claim of the Lemma follows after one makes a convenient choice of a = 3, 
b = 1/2 + 1/a. For this choice, A is close to its minimum for small values of p — 2. □ 

We recap our procedure upto this point: one first finds a ball in L P (C) around 
the approximate solution h* which is mapped into itself by the operator T v \ the exact 
solution g u+v of the Beltrami equation v = {y + f])9z +v satisfies the bound from the 
above Lemma. 

8. Estimates on the solution of the Beltrami equation in C 

Our next step is to estimate the error introduced in ignoring the part of the differential 
in a small neighborhood of infinity. As explained in the Introduction, define 



The solution of the Beltrami equation on a compact is a homeomorphism, hence 
the continuous inverse of g u+v is well defined. The differential (3 is again supported 
in a neighborhood of infinity. Here, rather then getting estimates on the size of this 
neighborhood from general considerations, we can use our knowledge of g", together 
with estimate ()7.1|) . to obtain a bound on r such that g u+r] (B>ji) contains D r . This 
amounts to evaluating on any collection of standard sets (disks) Zk in C that cover 
the circle {z G C : \z\ = R} and choosing the infimum of \g"(Zk)\. 

It follows immediately that the differential (3(z) = (3{R 2 / z)z 2 / z 2 where (3 is defined 
in ()8.1)1 is compactly supported in a disk of radius R 2 /r. In the following Lemma and 
in the rest of the paper, IHI^ will signify the L p -norm over a measurable set E. 

Lemma 8. Letj G Loo(C) be compactly supported in C\©r ; and satisfy ||7||oo < K < 1. 
Furthermore, let g 13 be the solution of the Beltrami equation g§ = j3gP where /3, given by 
(jff.ijl and j\4-^ , is compactly supported in H) R 2/ r . Then for all p > 2 such that KC P < 1 




-i 



(8.1) 



9 




z\ < C\z 



l-2/p 



where 



C = n 



1/p AK 



R 4/p(2 - KC P ) 
r 2 /p(l - KC P ) ' 



(8.2) 



Proof. Notice, that 



11/311 



oo 



oo 




oo 




oo — 



< K. 
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By a standard result from the theory of quasiconformal mappings (see, for example 
(Carleson and Gamelin 1991) ), if h is the fixed point of the operator T§, then 



h 



< 



(ttR 



4\Vp 



-Vp(1-KC p 



for all p for which KC P < 1 . By an argument similar to that of Lemma [7| 



I/O) - z\ < A 






h + 1 


\z\ 






oo 




D R 2 /r ,P 



ll-2/p < ^l/p^tf 



i? 4 /p(2 - 
r 2 /p(l - KC P 



\Z\ v . 



□ 



We can now complete the proof of Theorem El 



Proof of Theorem^ Let (3 be as in (JHTTj) and let, as before, (3(z) = f3(R 2 / z)z 2 / z 2 . Then 
the following bound on the solution g@ of gf = (3g@ is immediate: 



\A*)-A 



R 2 



g^R 2 /z) 



- - Ar 2 /z) 



< c 



R 2 



l-2/p 



gP{R 2 lz) 
C\z 



l+2/p 



C 



#il-2/p R 4/p _ C\ Z \ 2 /P 



Now, we can demonstrate the claim of the Theorem: 

\gP°g u+v (z)-g:(z)\ < \g^(z)-g:(z)\ + \g^{z) - / o g»+*(z) \ 

<Ac'\z\^ I C \9 U+V (*)\ 1+2/P 

R a /p -C\g^(z)\ 2/p 

< ^u.i- a /P + Cmz)\ + Ae>\z\^] l+2/p 

R 4lp ~ C [\gl{z)\ + Ae'\z\ 1 - 2 /pf p ' 

Notice, that the first term is e'-small, while the second term is l/r 2//p - small. 



□ 



9. Completion of the fundamental domain 



In this Section we will show how one can obtain bounds on the essential supremum and 
the Lp-norm of the Beltrami differential (jH.Hj) in the neighborhoods of zero and infinity. 
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As we have already mentioned in Section El the parts of the boundaries of the 
fundamental crescent adjacent to the indifferent fixed point can be chosen as a piece of 
an internal ray and its preimage under the branch of f~ q " that fixes 0. We will now 
briefly describe the construction of these boundaries and estimate the absolute value of 
the Beltrami differential for small values of the parameter. 



Let / be an analytic map on C with a fixed point at such that A = f /'(0) = e 
where 9 is of bounded type. By the Siegel Linearization Theorem, there exists a unique 
conformal local change of coordinates <f, normalized so that ip(0) = and y?'(0) = 1, 
that maps the disk © s of radius s = |<£> -1 (1)| onto the Siegel disk and conjugates / to a 
rotation by A: f(<p(z)) = (p(Xz). By the classical Koebe Distortion Theorem 

\z/s\ 

(l-\zJW 

for all z G D s . Therefore, if z G © p , p < s, then 



\<p\z)-l\< 



\<p(z) - z\ < c\z\ 2 , \ip\z) - 1| < c\z\, (9.1) 

where c = l/s(l — |p/s|) 3 . 

Choose a sufficiently small p > and define g to be the largest parameter value 
satisfying \ n (g) G ip(dB> p ). Define d = v ?_1 (A n (f?)) and consider a ray passing through 
the point A n (f>): the image ip(dt) of the short line dt, t G [0, 1]. The preimage f~ q "((p(dt)) 
of the ray passes through the point / _9,I (A n (^)) - - it will serve as a completion of 
the other boundary curve of the fundamental domain. As before, it is convenient to 
parametrize this part of the boundary in terms of the radial coordinate r. Here, we 
will make the choice t = r u /g u , the positive power v will be specified later. As we have 
demonstrated, the coordinate (p is close to the identity on D p , therefore the curve <p 
defined by setting <p(r) = tp(dr u / g u ) satisfies the following set of equations for all r < g: 

(p(r) = — r"(l + 7(r)), where |7(r)| < c\d\, 
Q 

m'( r ) = v— r v -\l + /3(r)), where W(r)\ < c\d\, 
g v 

\~ Qn d 

f qn ((p(r)) = — - — r I/ (l + 7(r)), where |7(r)| < c\d\, 
Q 

(f~ qn o$)'(r) = v^^r v - 1 {l + 0(r)), where |/3(r)| < c\d\, 
Q 

where, as before, A = e 2md . 
Next, 
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air) = In 



Qn 



g u d(l + 7(r)) 
is close to In (—a qn /d) — b for r < g. Similarly, 



r^r^r))) = - 
la 



air) — z/ln 2iri(—q n 8 mod 1) 

Q 



where 



rr(r) = In ( A~ g "^j — 



- 6 



is again close to In (—a qn /d) — b. Finally, 

cr(r) — 2m(—q n 6 mod 1) 



v(-<l>)+7r 

ZTT 



+ 



i-<i<-*)-£ 



cr r 



ln^ 
/ I. 

ia 



for small r, and, therefore, its complex dilatation on the complement of the non-negative 
semi-axis is given by 



\x{re 



a 2i<t> 



rd r g n (r,(p) +id^g n {r, 0) 



where 



a(r,(f>) 



rd r g n (r, 0) - id^g n (r,<f>) 

c 2i^ u ~ mod X ) ~ a ( r ^) ~ 2^( a ( r ) ~^( r )) 

v + {-q n d mod 1) - a(r, 4>) + ^(cr(r) - 5"(r)) ' 

(^(-^) + ^) ™'(r) + (l " V{-<I>) ~ ^) rcr'(r). 



(9.2) 



(9.3) 



Notice, that if z/ is chosen to be (—q n 9 mod 1), then the absolute value of this 
dilatation is close to 0. To provide a bound on this dilatation for all r < g, one needs 
to know functions ra'(r) and rcr'(r), as formulas ()9.2|) and ()9.3|) indicate: 



rcr'(r) 



ra r 



d(l + 7(r)) 



d(rV^)(l + 7(r))-a g7 , 
d(l + 7 (r)) 



7/ 



^ , Sn(/^W -7W) 



v- 



d\-i-{r u / g v )(l + 7(r)) - a 9n 



£?" (i(l + 7(r)) 2 
A _9 ™i/ h ^ 



d(l + 7 (r)) s 



It is not hard to evaluate functions a, o and a on a standard set / D (0, g), and 
therefore, obtain a bound on the essential supremum of the Beltrami differential ()9.2|) 
for r < g. 

As a next step, we will describe the boundary curves of the fundamental domain 
for very large parameter values. Let z = a qn + w, £ = f% n (a qn ), and let f(w) = 
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f Qn (a qn + w) — a qn . By the Koenigs Linearization Theorem if |£| ^ 0, 1 then there exists a 
unique local holomorphic change of coordinates ip, given by ip~ l {w) = linifc MOO (, k f~ k (w), 
such that ip(0) = 0, tp'(0) = 1 and f(ip(w)) = ip(£z). The latter limit is uniform on any 
disk D s such that < c|iu| for all w G B> s and some fixed constant c satisfying 

g a <ier 1 <5- 

As before, fix a small p < s, then according to the Koebe Theorem, there exist a 
constant c = l/s(l — |p/s|) 3 such that 



(w) — w\ < c\w\ 2 , \ip'{w) — 1| < c\w\ 



(9.4) 



Define i? to be the smallest parameter value satisfying X n (R) £ </?(<90 p (a ?r J), 
ip(z) = ip(z — a qn ) + a qn . Define d = </? _1 (A n (.R)) and consider the logarithmic spiral 
(d - a qn )e vhl ( r / R )/( 27r ) + a qn passing through the point ip~ l (\ n (R)) when r = R. Here v 
is some complex number, to be specified later, with a negative real part. The image 



<p(r) 



^{(d-aje^*) 



of this spiral is a simple smooth curve contained in (p(D p (a qn ))) for all r > R. For the 
sake of notational compactness below, we will use the same notation as in the case of 
the indifferent fixed point: of course, the functions and constants appearing in the two 
cases are completely unrelated. 



ip(r) = (d-a qn )e^ n R(l + -f(r)) + a qn , where \^(r)\ < c^d - a, 
(p\r) = (d-a qn )^e^ hl i(l + (3(r)), where 



Qn I 1 



Ml < c 2 \d-a qn \, 



and 



d — a n 



e 2* ln k(l + 7(r)) + a qn , where |7(r)| < d 



(r 1 o^)'(r) = ^ 



d — a qn e~- z " " 



In 



Next, 



_-i 



Mr)) 



-\r qn {^{r)) 



1 

ia 
1 

ia 



2nr 



In 1 - 



(l + /3(r)), where 



< c 2 "- 



\d — a q 
d~a qn 



(d-a qn )e^ la n(l + ^(r)) + a qn 



-b 



In 1 - 



r 1 (^-sJ^ ln ^(l + 7(r)) + S 



therefore 



<? n (r, 0) = r- 1 (#r)) + 1 ^(-0) + |-) ((7 ( r ) _ InO, 
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where 



a(r) = In 



l + 7(r) (d-a qn )(l + ^(r))e^ ln R +a, 



(9.5) 



1 + 7 ( r ) £-i( d _ + 7 ( r )) e £m£ + flg 

The complex dilatation of g n outside of the non-negative semi-axis is given by 

dre*) = e 2 ^ — ln £ +27rr " a ( r ) +r ( 27rr ?(-^) + ^) (a(r) - a(r)) +io-(r) g ^ 
v+i ln^ + 27rra(r) +r (27r?7(— 0) + 0) (a(r) — a(r)) — ia(r) ' 

where a and a, the r-derivatives of ior _1 (^(r)) and iar^ 1 (f~ qn (y5(r))) , respectively, are 
given by 

- f v l\r) (rf-a g Je^ ln i(f(r) + (l+7(r))^) 

" lrJ l + 7(r) (rf- S J(l+7(r))e^ ln i + a gn ' K ' } 



af _, r 1 (^-sJe^ ln i(Y(r) + (l + 7(r))^) 

l + 7(r) e- 1 (^-«,J(l + 7(r))e^ ln i+a ? 



Formulas ()9.5|) - (|9.8JI . allow one to estimate supremum of the Beltrami differential 
()9.6j) for r > R. We would like to emphasize that the absolute values of functions a, a 



and a are small for all r > R if d — a qn is small. 



The constant v can be now conveniently chosen to be equal to «ln£. 



10. Bounds on the Beltrami differential 

Finally, we would like to provide an estimate on the supremum and the L p norm of the 
error of the Beltrami differential for intermediate parameter values: r e [g, R] . 
Write the Beltrami differential (|3.3J) as 

aire 1 *) = e 2 ^ (r) + CM) + C(r) (10 1) 



where 



^ [Dt- 1 (r fe (A n (r))) Df~ qn (A„(r)) - Dr" 1 (A n (r 



MH = r 

27T 

«r) = ^r-\r qn {Ur))) ~ ^(A^r)), 

C(r, 0) = »7(-^)C-(r) + (1 - r7(-0))C+(r), 

C-(r) = r\' n (r)DT-\r qn {Ur)))Df- qn (MO) 
C+(r) = rAUr) J Dr- 1 (A n (r)) . 
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2tt 



MO = ;f- / Kre ta )e- tka da = — 



1 



2tt J_ v n(r)a + ((r, a) - £(r) 



e^ a da. 



Therefore, one obtains after integration: 



H2{r) = 1 + 



7r/t(r) 



In 1 + 



7T 



a+(r) 



In 1 



7T 



«_(r) 



/i 2 _ fc (r) = r ^ [Ei (ifc (tt + a+(r))) - Ei (ifca+(r))] 

+ ~ Urj [Ei (ifca_(r)) - Ei (ik (a_(r) - tt))] 

7T/c(r) 

= f-l) fe ^^ [ExpEi (ifc (tt + a+(r))) - ExpEi (life (a_(r) - tt))] + 



+ 



e(0 

7r/c(r) 



7r/c(r) 

[ExpEi (ifca_(r)) - ExpEi (i/ca + (r))] , 



dcf 



where ExpEi (z) = e z Ei(z), function Ei being the exponential integral, and 

a (r) - C ~ (r) -^ (r) a (r) - C+(r) = ^ 
a - (r) " «(r) ' + () ~ «(r) • 

A simple inspection shows that ai+(r) + 7r = a_(r) — n, therefore 



1 + (a_(r))- In (a + (r))], fc = 0, 

^ [ExpEi (jfea_ (r)) -ExpEi (ifca+(r))], fc ^ ' 
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(10.2) 



Function ExpEi is holomorphic in the complex plane with a branch cut along the 
positive real axis. Of special interest to us is the fact is that the exponential integral 
has the following representation at zero: 



Ei(z)=j + ]n{-z)+J2?- + E n (z), \E n (z)\ < 



i=i 



(n + l)!(ra + l)' 



10.3) 



where 7 is the Euler constant. Furthermore, function ExpEi itself admits the following 
asymptotic series at infinity (cf (Lebedev 1972)): 



ExpEi (z) 



.i=0 



(10.4) 



Renormalization for Siegel disks and Measurable Riemann Mapping Theorem 



23 



where 



\R n (z)\ < [ V (Rez)(\z\ - \lmz\) + 1] 



(10.5) 



We have used series (jl(J.3j) and (|1U.4|) to provide bounds on the function ExpEi , 
and, therefore, bounds on the coefficients fi2-k, k ^ 0, of the Beltrami differential. 
Moreover, since the typical values of a± are large, asymptotic series (jll).4J) allows one to 
compute a bound on the supremum of the Beltrami differential on \ D e . To establish 
such bound, first, observe that 



TCK(r 



7rk 



l ~ kt^~ pj^ + ]J%r + -j^r + R 5 (ika + 
ika + {r) 

3i 



(10.6) 



for all fc ^ 0. Therefore, /i 2 _ fc + fi 2+ k = A{r)/k 2 + B(r)/k A + (iT(r) + R + {r))/k 6 and 
H2-k ~ = iC{r)/k + iD{r)/k 3 + iE{r)/k 5 + {R~{r) - R + {r))/k 6 , where 



A(r 
B(r 
C{r 
D(r 
E(r 



2£(r) 



1 



7T V(C+(r) -£(r))a+(r) (C-(r) - £(r))a_(r) 



6£(r) 



* V(C-(r) -eW)ai(r) (C+(r) - £(r)K(r) 



2£(r) 



vr VC + (r)-^(r) C-(r) - £(r) 



4£(r) 



* V(C-(r) -C(r))a 2 - (r) (C+(0 " 



8£(r) 



and 



V(C+W -{(r))4(r) (C-(r) -eW)ai(r) 
f(r) / i? 5 (^'Q; + ) i? 5 (^o;_) 
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Next, 



\fi(z)\ 



< 



< 



fc=i 



Ata(r) + y](At2-fc(r) + H2+k(r)) cos A;0 - i((j, 2 -k(r) - H2+k(r)) sin /c0 
fc=i 

/ \ v-^ , cos kd) „. . cos kd) .smkd) . sinfc0 

/x 2 (r) + ^ A(r)^^ + ^(r)^^ + C(r)— r-^ + £>(r)^r^+ 



fc=i 



+£(r) 



A; 2 



A; 4 



fc 



A; 3 



—ik<f> 



D ik<f> 



sinkrfj _ , 

— — - + i2 r — — + i? + r — 

A; 5 A; 6 A; 6 



;io.7) 



Therefore, sup 2e Dfi ^ Dg | is bounded by the supremum of the following expression 
on [g,R] x (0,2tt): 



/j 2 (r) + A(r) 
+ D(r) 



7T 2 7T0 2 

y ~ t + t 



7T 



7T0 2 3 



6 



12 



+ 5(r) 
+ £(r) 



7T 4 7T0 2 7T0 



+ 



90 12 12 48 



7rV 7T 2 (/> 3 7T</> 5 4 



+ C7(r 



7T — 



90 36 



+ 



48 240 



7T 



+ 945 (l ^ (r)l + |i?+(r)l) - 



(10.8) 



Lastly, we will estimate the accuracy of the approximation of the Beltrami 
differential with a finite Fourier series. Fix some sufficiently large even natural M 
and represent /x on U> R as v + 77, where v is supported in \ U> g and r\ = ji on D e . 
Furthermore, we can assume that v is representable as a finite Fourier series on D fl \D e , 
and r] is the "higher order error" on D R \ D e : 

M+2 



k=2-M 
oo 



l-M 



V[ re^)= 



e ik<t, + n k (r)e ik<t > = 771 (re**) + r? 2 (re^), 



77! (re**) =e 210 ^ 

fc=M+l 



fc=M+3 fc=-oo 
oo 

2i</> 



A(r) S(r) 
~fc 2_ + _ fc 4_ 

ik(f> 



COS /C0 + 

,ik(f> 



C(r) _D(r) £(r) 



A; 



A; 3 A; 5 



sin A;0, 



?7 2 (re^)=e^ £ iT(r)^ + i? + (r) 



fc=M+l 



A; 6 
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To obtain a bound on the L p -norm of 77 on IB># \ D e , one can use the fact that 
v\\mAn:.v < su P^d^ W^I^IMI^Wo for p > 2, together with the following 



bound on the squares of the L 2 norms of r]x and r\2 '■ 

„^> oo 1— M 



- k=M+3 k=-oo 

/R 00 00 
|/i 2+fc (r)| 2 + Y \^2-k{r)\ 2 rdr 
^ k=M+l k=M+l 
rR 00 

= 71 / l^ 2 - fc ( r ) + + \l*2-k(r) - ^2+k{r)\ 2 rdr 

J p i — 1 1 



'e k=M+i 

N-l 



- 0^! 2 Cm(2) + [|^| 2 + 2Re (QA)] Cm(4) + 



i=0 

|2 



2 

[| A| 2 + 2Re (AA) + 2Re (A.A)] Cm(6)- 
[| A| 2 + 2Re (AA)] Ca/(8) + | A| 2 Cm(10)] (p 2 - p 2 +1 ), (10.9) 

\m\\l<2n f \dr [\R-(r)\ 2 + \R + (r)\ 2 ]( M (12) 



N-l 



= ttCm(12) J> 2 ~ pIi) [m 2 + \Rt\ 2 } , (io.io) 

i=0 

where Cm(^) are the (M + l)-remainders of the Riemann zeta function: Chii 71 ) = 
Y^=M+i^~ n -> an< l ^i, Bi, Ci, A) A and i?^ 1 are bounds on the corresponding functions 
on intervals (p^, Pi+i) that cover (g,R), Po = Q (see Appendix for details). We would 
like to emphasize that the remainder of the zeta function Cm (2) decreases very slowly 
with a growing M. The appearance of this function in the leading term of the norm of 
r) necessitates working with high order truncations of the Fourier series. 

Naturally, bounds (|l().8jl . ()10.9|) and ()10.10|) can be implemented on a computer. 

The Lp-norm of rj on 3 e can be bounded by sup 2gE)p \p.(z) \ (vrf? 2 ) 1 ^, where 
sup 2gD is computable with the help of formula ()9.2|) . 



11. Appendix. Standard sets 

Over the course of the last two decades computer-assisted proofs have become a standard 
mathematical tool in dynamics since they first appeared in O. E. Lanford's (1982, 1984) 
pioneering works in the early 80-s (also, cf (Eckmann et al 1982), (Eckmann et al 1984), 
(Mestel 1989), (de la Llave and Rana 1990), (Stirnemann 1994) and (Koch 2004a) for 
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some other important computer-assisted results). Here we will give only a very brief 
description of some ideas: An interested reader is referred to an excellent review (Koch 
et al 1996) on the subject. 

Given a set X, denote by V(X) its power set. Let / be a map from a subset 
D f C X to another set Y. We say that F : D F C V(X) ^ V(Y) is a bound on / if 
f(z) G F(Z) whenever z G Z G Dp. It is customary to construct such bounds within 
the class of maps F : std (X) i— > std (Y), where std (X), the collection of standard sets 
in X, is a conveniently chosen subset of V(X). 

11.1. Interval arithmetics in C 

A computer implementation of an arithmetic operation r^r-i (# is +, — , * or /) on 
two real numbers does not generally yield an exact result. The "computer" result is 
a number representable in a standard IEEE floating point format (cf (IEEE 1985)). 
Such numbers are commonly referred to as "representable". However, it is relatively 
straightforward to find the smallest representable number larger than ri#r 2 and the 
largest representable number smaller than ri#r 2 , if such numbers exist. If either of 
these numbers, denoted r> and r<, does not exist then the pair (r 1? r 2 ) is considered to 
be outside of the domain of the map (ri,r 2 ) i— > ri^r 2 and the program is aborted. If 
such bounds do exist, then one can take the interval /[ bound on the result 

of the operation. Thus, one is naturally brought to the idea of constructing bounds on 
functions on real numbers as functions on the set of all representable intervals: 

std(R) = {I[x,y} G R : x, y - representable}. (H-l) 

To obtain a bound on a transcendental function (exp, cos, sin, In and others) one 
can use their Taylor series together with a bound on the remainder. Some algebraic 
functions, like roots, can be bounded with the help of these transcendental functions. 

Arithmetic operations on complex numbers are reduced to those on reals, and so is 
the evaluation of bounds of functions on C on standard sets in C. We will choose these 
standard sets to be open disks in the complex plain with representable centers and radii: 



std (C) = {B[z, r] = U) r (z) C C : z = x + iy; x, y, r — representable}. (H-2) 
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11.2. Operations on standard sets for analytic functions 

Given a function / G C p we will bound it by a standard set of the following form: 
F (B, I g , I h ) = If G C p : f(z) = p{z) + zg{z) + z N+1 h(z); 

fc=i j 

where A^ is a sufficiently large fixed integer, B is the shorthand for the array of 
Bk G std(C), k = 1...N, and I g , Ih are in std(lR). The function g(z) is commonly 
referred to as the "general error" (cf (Koch et al 1996) for the explanation of the necessity 
of this term), the function h(z) is usually called the "higher order error". 

Standard sets (jll.3j) are well suited for obtaining bounds on addition, 
multiplication, composition and differentiation of functions in C p . As an example, 
we will consider bounds on multiplication of two analytic functions. Let f\ and /2 be 
in C p . Given two bounds on these functions, F x 3 /i and F 2 3 f 2 , we would like to 
obtain a standard set F G std(C p ) which contains fx ■ f 2 - Let pi = J2k=i a k zk an d 
p 2 = J2k=i bkZ k - Let p\{z) ■ p 2 {z) = p{z) + z N+1 r(z) where the order of p is no more 
than N. Then 

(fi ■ f2){z) = p(z) + z \px(z) ■ g 2 (z) + g x (z) ■ p 2 (z) + z ■ g x (z) ■ g 2 (z)} 

+ z N+1 (r(z) + z N+l ■ h(z) ■ h 2 (z) + ( Pl (z) + z 9l (z)) ■ h 2 {z) 

+ h(z) ■ (p 2 (z) + zg 2 (z))) . (11.4) 

Here, the terms in the brackets constitute a general error, while those in the 
parentheses are the higher order error terms. One knows the bounds on all objects 
that enter these error terms, therefore one can find the interval I g and Ih that contain 
the norms of the general and higher order terms of fi ■ f 2 . 

We would like to note that although our programs include all the necessary routines 
that supply bounds on operations with standard sets (jll.3j) . these routines have not been 
fully used in the proof of the Main Example. This Example serving only demonstrative 
purposes, the analytic function appearing in it is a quadratic polynomial. Operations 
with this polynomial do not require estimates on general and higher order terms (if 
is chosen to be sufficiently high). 

11.3. Standard sets and operators in L P (C) 

We have adopted the strategy of approximating a function hire 1 ^) = Y2T=-oo hk(r)e ikt ^ in 
L P (C) by a finite Fourier series with piecewise-constant compactly supported coefficients. 
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Given a collection of real numbers = r < r\ < . . . < r^-i < < r^+i < . . . < tm 
that specify a radial grid, we will denote the midpoint of the interval [r m , r m +i] by g m ; 
Q-i will be by definition. We will make the following choice of the standard sets in 
L P (C) associated with a grid {r m }m=o- 

{M 2 
h = s + ee L p (C) : s(re^) = ^ s k {r)e ik(t >, s k (r) e S Km for r G 
fc=-Mi 

[gm-i, £m),ra < and Sfc(r) = for r > ||e|| p e 7 e | ■ (H-5) 

Here, Mi and M2 are fixed positive integers, and S is the shorthand for the array of 
S Km e std (C). 

Let h = s + e e L P (C) lie in a standard set H (S,I e ). We will now describe the 
standard set that contain T[h] and With this goal in mind, we will describe the 

action of the Hilbert and Cauchy operators on piecewise-constant functions. 

Let s be a compactly supported function in L p (C) given by a finite Fourier series 
with piecewise-constant coefficients: 



) = £ ..(.V* .,(r) = { Sk "" ' £ [<? ™-;' £ - N ' (11.6) 



fc=-Mi 



The action of the Hilbert operator on such function again can be represented as a 
Fourier series: 

M 2 

T[ a ](re*)= £ c k (r)e ik *. (11.7) 

fc=-Mi 

The coefficients of this series are given by the following set of expressions (cf 
(Gaydashev and Khmelev 2005) and (Daripa and Mashat 1998)): 

1) r = 

' c (r 1 )-s 2 (r 1 )-2s 2 (r 1 )\nf o , k = 0, 
0, k ^ 0. 

2) r-E [Qm-uQm), 0<m<N, 



Cfc(0) 



Cfc(r) 



(^Tt) ( C fe( r ™+l) ~~ Sfe+2(^m+l)) - 4,m(r) + Sfc +2 (r), fc > 1, 

c (r m+ i) - s 2 (r m+1 ) - 2s 2 (r m+1 ) In ^ - 2s 2 (r m ) In ^ + s 2 (r), fc = 0, 

(^wT) ( C fc( r ™~l) - S fe+2(^m-l)) + 6fc,m(»") + S fc+2 (r), < 0, 
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dk, m {r) = 2 
bk, m (r) = 2 



fc + 1 
fc 

fc + 1 
fc 



Sfc+2(^ m+l} 
•Sfc+2(j m-1, 



cm m+1 



TO— 1 cm— 1 



3) r6 [ftn-i, g m ), N <m< M 

Cfe(r) 



■Cfc(r ri 



3 HfcV' m-i;> 

'm-l 



+ s k+2 {r n 

fc^O, 
k < 0. 



i- + 

£m 
c*m— 1 



To find the standard set H(C, I[0,0]) that contains T[s] one has to identify the 
standard sets C^.m such that Cfc(r) G Cfc jTO whenever r G [f? m -i,f?m)- This amounts to 
evaluation of Cfc(r) on an interval I m G std (R) containing [g m -i, g m ). 

Finally, the interval lT[ e ] that contains ||T[e]|| p can be readily found with the help 
the Calderon-Zygmund lemma: It is the unique standard set with the smallest length 
in std (R) that contains all real numbers C p x for all x G I e . 

In a similar fashion the action of the Cauchy transform on a piecewise-constant 
function (jll.6j) is given by the following set of equations: 



1) r = 

Cfc(0) = 

2) r G [g m -i,g m ), 0<m<N 



c (ri) - 2s 1 (r 1 )g - 2s 1 (r 1 )(r m - g ), k = 0, 
0, fc^O. 



Cfc(r) 



V fc 

^7 J c fc (r m+1 ) - d htm (r), 



ci(r m +i) - 2rs 2 (r m ) In ^ - 2rs 2 (r m+ i) In 



-) c fe (r m _i) + 6 fe , m (r), 



fc > 1, 
fc = 1, 

fc <= 0, 



where 



d k , m {r) 



h,m(r) 



2r 



1 - fc 
2r 



Sfc+i(r r , 



em 



m+1 



1 - fc 

3) r G [ft n -i, ft»), N <m< M 

Cfc(r) 



„i fc ^,1 fc 

cm— 1 ' m— 1 



,1-fc 



,1-fc 



fc^O, 



cm 
cVt-1 



o, 

-c fc (rm-i), fc < 0. 
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The bound on P[s] is found similarly to the analogous estimate for the Hilbert 
transform. Notice, that we do not need the bound on P[e\. According to Theorem El 
one only requires bounds on the action of the Cauchy transform on the approximate 
fixed point h* of the iteration scheme. Such approximate fixed point is always given as 
a finite Fourier series. 
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